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We present a systematic study of the frequency band structure of acoustic waves in crystals 
consisting of nonoverlapping solid spheres in a fluid. We consider colloidal crystals consisting of 
polystyrene spheres in water, and an opal consisting of close-packed silica spheres in air. The opal 
exhibits an omnidirectional frequency gap of considerable width; the colloidal crystals do not. The 
physical origin of the bands are discussed for each case in some detail. We present also results on 
the transmittance of finite slabs of the above crystals. 



PACS numbers: 43.20.- 
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I. INTRODUCTION 

Colloidal suspensions of monodisperse spherical parti- 
cles, with a diameter betwreen 1 nm and 1 /im, in liquids 
or gases self assemble onto three-dimensional simple lat- 
tices whose lattice parameter can be easily tailored, pro- 
viding new opportunities for fundamental as well as for 
applied research (see, e.g., Ref. |l| and references therein). 
The optical properties of these colloidal crystals are be- 
ing investigated by many research groups because they 
seem to be good photonic band-gap materiajsjJjtut they 
have interesting acoustic properties as welLHuaB More- 
over, experiments relating to the reflection, transmission 
and absorption of ultrasonic waves by colloidal crystals 
can be very useful for the characterization of such sys- 
tems, provided the means exist for a proper theoretical 
analysis of the experimental data. Finally, theoretical 
results for colloidal crystals offer a starting point for the 
understanding, at a semiquantitative level, of the prop- 
agation of acoustic waves in corresponding random me- 
dia c^^t^ting of monodispersed spherical particles in a 

The vibration modes (normal modes of the elastic field) 
of a phononic crystal, by which we mean a composite ma- 
terial whose density p(r) and Lame coefficients A(r) and 
//(r) vary periodically in space, are Bloch waves with a 
corresponding frequency band structure which is analo- 
gous to that of electrons in ordinary cisHtals--apd|-e|lcctro- 
magnetic waves in photonic crystals .B'BBclMO'tHl With 
an appropriate choice of the parameters involved one may 
obtain phononic crystals with absolute (omnidirectional) 
frequency gaps (phononic gaps) in selected regions of fre- 
quency. An elastic wave, whose frequency lies within an 
absolute frequency gap of a phononic crystal, incident 
on a slab of the crystal of certain thickness will be re- 
fiected by it, the slab operating as a perfect nonabsorb- 
ing HHpmr of elastic waves in the frequency region of the 
gapJldllil It may be possible, for example, to construct 




FIG. 1: Scattering of a plane acoustic wave by a polystyrene 
sphere in water: total scattering cross section. 



in this way vibration-free cavities which might be very 
useful in high-precision mechanical systems operating in 
a given frequency range. 

In this paper we investigate in detail the acoustic prop- 
erties of fee colloidal crystals of polystyrene spheres in 
water and of close-packed silica spheres in air (opals). 
Given that these systems have a characteristic length 
scale of the order of 1 /im, they should exhibit interest- 
ing acoustic properties at ultrasonic frequencies of a few 
hundred MHz to a few GHz. It should be noted, however, 
that our results apply to different regions of frequency of 
the acoustic field provided that the size of the spheres and 
of the unit cell are scaled accordingly, and provided the 
elastic coefficients can be taken as constants independent 
of frequency over the said regions. The calculations of the 
frequency band structure and of the transmission coeffi- 
cient of acoustic waves through a slab of the material 
were done using the layer-multiple-pscattering formalism 
we have developed for this purposcHil A formalism-aJong 
the same lines has also been published by Liu et al. .113 We 
have already demonstrated the efficiency of this method 
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FIG. 2: (a) The phononic frequency band structure normal 
to the (111) surface of a fee crystal of polystyrene spheres 
in water with / = 30%. The solid/dotted/dashed lines refer 
to bands of A1/A2/A3 symmetry, respectively, (b) Transmit- 
tance of an acoustic wave incident normally on a slab of the 
above crystal consisting of 32 planes of spheres parallel to the 
(111) surface. The positions of the first two resonant modes 
of a single polystyrene sphere in water are indicated in the 
margin. 

in relation to solid-solid composites. The present paper 
shows that the method applies equally well to phononic 
crystals consisting of nonoverlapping solid spheres in a 
fluid host; it appears that the plane- wave method for 
calculating the frequency hand structure of such systems 
has convergence problems.EI Besides, we shall pay particu- 
lar attention to the physical origin of the different modes 
of the acoustic field in the systems under consideration. 



II. POLYSTYRENE SPHERES IN WATER 

We consider a model colloidal crystal consisting of 
polystyrene spheres in water. The mass density and the 
sound velocities for polystyrene are: pp = 1050 kg/m"^, 
cip = 2400 m/s, Ctp = 1150 m/s. For water we have: 
Pu) — 1000 kg/m^, = 1480 m/s. We begin by consid- 
ering the scattering of a harmonic plane acoustic wave 
(a longitudinal wave) of angular frequency w, by a sin- 
gle polystyrene sphere of radius S in water; the water 
extends over all space. The incident plane wave can 
be written as a sum of spherical waves associated with 
the spherical harmonics Yjm, where I = 0, 1,2, ••• and 
m = —I, • • • , 0, • • • , Z, as usual.llil A spherical wave of 
given (l,m) scatters independently of the others, of dif- 
ferent {l,m), because of the spherical symmetry of the 
scatterer; therefore the total scattering cross section is 
the sum of partial (Z, m) cross sections, with Z up to a 
maximum Imax depending on the size parameter, ojS/cw, 
of the sphere. We assume that waves with / > Imax do not 
scatter from the sphere and do not contribute to the total 
scattering cross section. Of course partial cross sections 



FIG. 3: The phononic frequency band structure normal to 
the (111) surface of a fee crystal of polystyrene spheres in 
water with / = 30% (only bands of Ai symmetry are shown) , 
calculated (a): with Imax = 1 and (b): with Imax ~ 6. 



of different m (but of the same /) are equal because of the 
spherical symmetry of the system. In Fig. |l] we show the 
total scattering cross section for a plane wave scattered 
by a polystyrene sphere in water as a function of LuS/cm. 
The peaks correspond to resonant modes of the acous- 
tic field about the sphere with an angular distribution of 
the displacement field at the surface of the sphere deter- 
mined by the spherical waves which contribute to these 
modes. We find that the resonance at loS/cw — 1.49 
is an Z = 2 resonance: the displacement field associated 
with it is made up almost entirely (98%) from I — 2 
spherical waves; the one at loS/cw = 2.19 is an / = 3 
resonance: the corresponding displacement field is made 
up mostly (90%) from I — 3 spherical waves; finally the 
field associated with the resonance at uS/c^ — 2.83 is 
made up mostly from I = 4 waves (65%) and / — 1 waves 
(22%). And we remember that there will be {21 + 1) res- 
onant modes of the displacement field corresponding to 
TO = — • • • , 0, • • • , Z of the same frequency ui, i.e. coi is 
{21 + l)-degenerate. 

We now consider a fee crystal, with lattice constant a, 
of nonoverlapping polystyrene spheres in water. We view 
the crystal as a sequence of (111) planes of spheres. The 
spheres of a plane are arranged on a hexagonal lattice 
of lattice constant ao — a^/2/2, defined by the primitive 
vectors ai — ao(l,0) and a2 — ao(l/2, \/3/2), in the xy 
plane (plane of the spheres). The {n + l)th plane along 
the z axis is obtained from the nth plane by a primitive 
translation, as, of the crystal. The planes are separated 
along the z direction by a distance d = a^/3/3. 

Fig. 11(a) shows the frequency band structure of acous- 
tic waves in this crystal normal to the (111) surface, ob- 
tained 
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when the fractional 
= 30%. In this case. 



by the method of Ref. 
volume occupied by the spheres is / 
the component of the reduced wavevector within the sur- 
face Brillouin zone (SBZ) of the fee (111) surface, k||, 



3 



equals zero. The symmetry of the bands ^^((jj; k|| = 0) is 
determined by the symmetry group for this direction (the 
Cat, group) and the corresponding propagating modes 
(Bloch waves) of the acoustic field belong to the dif- 
ferent irreducible representations of this group, denoted 
by Ai, A2 and A3 (see, e.g., Ref. |l6|); Ai and A2 are 
one-dimensional: the corresponding frequency bands are 
non-degenerate, and A3 is two-dimensional: the corre- 
sponding frequency bands are doubly degenerate. At the 
r point (k|| = 0, fcz = 0) the modes of the acoustic 
field belong to the irreducible representations of the Oh 
group associated with this point, as shown in Fig. ||(a). 
We see that the frequency bands of Fig. ^(a) appear in 



pairs: kz{u}]'k\\ = 0), —kz{i^\'k\\ = 0). However, this 
symmetry property is not valid for an arbitrary value of 
k|| , because the crystal under consideration, described by 
the Oh group, does not have a plane of mirror symme- 
try parallel to the (111) surface which would transform 
{kx, ky, kz) to {kx,ky,—kz). The above symmetry prop- 
erty holds only for k|| on the x axis [k|| — {k^, 0)], because 
a rotation through tt about the x axis is a symmetry op- 
eration of Oh and it transforms (fca;, 0, kz) to (fc^;, 0, — fc^). 

The (longitudinal) wavefield, of given lo and k|| , in the 
host region between the nth and the (n -|- l)th planes of 
spheres, can be expanded into plane waves propagating 
(or decaying) to the left and to the right as follows 



u(^^;k||) = V<^ u+ exp 



zK+ .(r-A„) 



u exp 



«K„, • (r- A„) 



(1) 



r 



with 



K± = (k|| + g, ± [(C^/C^)2 - (k|| + g)2] 



1/2 



(2) 



where g are the two-dimensional reciprocal vectors cor- 
responding to the lattice defined by ai and a2 above, and 
A„ is a point between the nth and (n -I- l)th planes. A 
generalized Bloch wave satisfies the equation 



ur 



exp (ik • a3) u^^ 



(3) 



where a3 = A„+i —A„ and k = (k|| , ^^(ci;; ky)) . We note 
that, although both longitudinal and transverse modes of 
the elastic field are considered within the (solid) spheres, 
only longitudinal modes exist in the (fluid) host region 
in the absence of viscosity. Therefore, in a binary com- 
posite of nonoverlapping solid spheres in a fluid host, 
where the solid component does not form a continuous 
network, there cannot be propagating transverse waves. 
Consequently, at the F point (k|| = 0, k^ = OHhe disper- 
sion curves [Fig. ||(a) and also Figs. |[ ^, and 0(a) below], 
show only one (longitudinal) branch starting from zero 
frequency, instead of the three (corresponding to both 
longitudinal and|-jt«Mp4eiSje-|Bipdes) appearing in solid- 
solid composites.BiByJl3y.y 

One can easily show that when k|| — 0, the g = 
component of the wavefield, described by Eq. (^, van- 
ishes (Ug^o = 0) for the modes of the A2 or the A3 
symmetry; only the Ai symmetry allows a non-vanishing 
g = component of the wavefield. Now assume that we 
have a slab of the crystal of finite thickness, i.e. N fee 
(111) planes of polystyrene spheres; between the spheres 
and to the left and right of the slab, extending to infin- 
ity, we have water. There, and for the frequency range 
which interests us here, the acoustic wavefield has only 
the one component g = 0. And therefore the external 
field couples with the field inside the slab essentially only 



through the g = component of the latter. Therefore 
an acoustic wave of given frequency, incident normally on 
the (111) slab of the crystal, will excite essentially only 
Ai modes of the crystal. If such a band does not exist 
at the given frequency the wave will be totally reflected 
by the slab. There will be no transmitted wave. The 
reader will see that this is indeed the case in the present 
instance: we show in Fig. |^(b), opposite the frequency 
band structure, the transmission coefficient of an acous- 
tic plane wave incident normally on a slab of the crystal 
consisting of 32 planes of polystyrene spheres in water. 
The transmission coefficient opposite the Ai bands ex- 
hibits the well-known Fabry-Perot-like oscillations due to 
multiple scattering between the surfaces of the slab; else- 
where it practically vanishes, and this includes regions of 
frequency where only bands of the A2 and/or A3 symme- 
try exist. Of course the above argument holds for normal 
incidence; off this direction the symmetry argument does 
not apply, and there will be some transmission if there 
are any propagating modes of the acoustic field at the 
given frequency (see below. Fig. |4|). 

From a practical point of view, once the frequency 
band structure and, if required, the transmission coef- 
ficient have been calculated, one has all that is necessary 
for a comparison with relevant experimental data. How- 
ever, it is worthwhile to look at the physics behind the 
band structure of Fig. |[ We note that the five narrow 
bands about toS/cw ~ 1.5 (we remember that the A3 
band is doubly degenerate) derive from the I = 2 reso- 
nances on the individual spheres, which interact weakly 
between them. The five-fold degeneracy of the reso- 
nance of the single sphere is split by this interaction 
in accordance with the lower (cubic) symmetry of the 
crystal field. Similarly the seven narrow bands about 
LuS/cw « 2.2 derive from the I = 3 resonances on the indi- 
vidual spheres. We observe the small hybridization gaps 
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FIG. 4: (a) The phononic frequency band structure associated 
with the (111) surface of a fee crystal of polystyrene spheres in 
water with / = 30%, for ky = |^(0.3, 0.1). (b) Transmittance 
of an acoustic wave incident with ky = |^(0.3,0.1) on a slab 
of 32 fee (111) planes of polystyrene spheres in water. 

opening up about LoSjCw « 1.5 and about LoSjc^ « 2.2, 
between bands of the same symmetry. 

There is also a mode of propagation of acoustic waves 
corresponding to almost free propagation in an effective 
homogeneous medium. In the absence of the resonances 
associated with the spheres, this mode of propagation, 
which has the Ai symmetry, would dominate the fre- 
quency band structure of the acoustic field as shown in 
Fig. 11(a). It was calculated by suppressing the I > 1 
resonances, i.e. cutting off the I > 1 spherical waves, 
which give rise to these resonances, from the spherical- 
wave expansion of the wavefield (see Ref. |Tl|). In the 
long-wavelength limit (u — > 0) we obtain a linear disper- 
sion curve, the slope of which gives an effective velocity 
of sound for the composite medium, c = 1589 m/s, which 
is in very good agreement with the result, c = 1566 m/s, 
of the effective-medium approximation]13 The small gap 
about LoS/cw ~ 1.5 in Fig. |^(a), is a Bragg gap; it is 
analogous to the small gaps one obtains at the edges of 
the Brillouin zone in the electronic band structure of free- 
electron-like metals. When the resonances on the spheres 
are allowed in, one obtains the band structure shown in 
Fig. |^(b) , with apparent hybridization gaps between the 
continuum band [of Fig. ^(a)] and resonance bands of the 
same symmetry (Ai). 

We believe that the disorder that is naturally there 
in a colloidal solution, does not eradicate the essential 
characteristics of the acoustic modes as calculated here, 
which are determined by the local environment about a 
sphere and less so by long range order. We would ex- 
pect the fine features associated with narrow bands to 
be smoothed out by disorder, but bands of modes sepa- 
rated by gaps could remain. This would be in accordance 
with the results of Brillouin-scattering experiments on 
disordered colloidal suspensions of monodispersed poly- 
methylmethacrylate (PMMA) spheres in transparent oil. 
It has been shown that in these colloidal suspensions dif- 
ferent longitudinal modes of propagation of the acous- 




FIG. 5: The phononic frequency band structure normal to 
the (111) surface of a fee crystal of polystyrene spheres 
in water with (a): / = 20% and (b): / = 50%. The 
solid/dotted/dashed lines refer to bands of A1/A2/A3 sym- 
metry, respectively. The positions of the first two resonant 
modes of a single polystyrene sphere in water are indicated in 
the margin. 



tic field exist which are separated by frequency gaps 
about tiicj-esonance frequencies of an individual PMMA 
sphere,clQcl which agrees, at a semiquantitative level, 
with the hybridization-induced gaps discussed above. 

In Fig. ^(a) we show an example of the band structure 
for a k|| 7^ 0, and next to it, in Fig. ^(b), the trans- 
mission coefficient of an acoustic wave, with the same 
kj[ , incident on a slab of the material consisting of 32 
fee (111) planes of polystyrene spheres. The thing to 
note is that all bands are active in this case, although 
the transmission coefficient is not always unity [compare 
with Fig. H(a)]. We remember that the incident wave 
cannot have a frequency smaller than Wi„/ = Cu, |k|| | . Fi- 
nally, it can be seen that now the frequency bands do not 
have the refiection symmetry found in the case of k|| = 
[Fig. ^(a)], for the reasons given above. 

We calculated the frequency band structure for a suf- 
ficient number of k|| points within the SBZ. There is a 
near-omnidirectional gap about luS/cw ~ 1.3; we tried to 
turn this gap into a proper omnidirectional gap by chang- 
ing some of the parameters, but we did not succeed. 

In Fig. H we look at the dependence of the frequency 
band structure on the fractional volume occupied by the 
spheres. We see in particular that the width of the 
resonance bands increases with /, apparently because 
the spatial overlap of the wavefield associated with reso- 
nances on neighboring spheres increases with /. 



III. OPALS 

In this section we consider a phononic crystal, an ar- 
tificial opal, consisting of two media (silica spheres in 
air) with density and velocity contrasts much higher than 
was the case in the crystal (polystyrene spheres in water) 
studied in the preceding section. Opalescent structures 



FIG. 6: Scattering of a plane acoustic wave by a silica sphere 
in air. Solid lines: partial scattering cross sections corre- 
sponding to I = 0,1,2. Dashed line: total scattering cross 
section. 



(usually referred to as synthetic or artificial opals) can 
be obtained from monodispersed silica colloids, e.g. by 
sedimentation in the gravitational field. In this way, fee 
arrays of silica microspheres in air, near the close-packing 
density (pj^ 74%), similar to a natural opal, have been 
obtainedjj Here we present some results on the acous- 
tic properties of such an opal. The mass density and 
the sound velocities for sihca are: ps = 2200 kg/m"^, 
cis = 5970 m/s, Cts = 3760 m/s. For air we have: 
Pa = 1.23 kg/m3, ca = 340 m/s. 

Proceeding as in Section ||, we first consider the scat- 
tering of a harmonic plane acoustic wave, of angular fre- 
quency u>, by a single silica sphere of radius S in air. 
In this case the total scattering cross section is a slowly 
increasing function of ivS/ca (see Fig. ||). However, the 
partial cross sections, corresponding to different values 
of /, exhibit well defined resonance peaks, of some width, 
as can be seen from Fig. ^ At uS/ca = 2.22 we obtain 
an I = resonance, at uS/ca — 3.88 an / = 1 resonance, 
and at LoS/ca = 5.25 an I — 2 resonance. 

We consider then an fee crystal of nonoverlapping sil- 
ica spheres in air with / = 74% and view it as a sequence 
of (111) planes of spheres. Fig. |^(a) shows the frequency 
band structure of the infinite crystal, for k|| = 0. The 
physical origin of the bands shown in Fig. 0(a), other 
than the fiat band about ujS/ca — 3.426, can be under- 
stood in much the same way as for the bands of Fig. ||(a). 
In the absence of hybridization we would have, extending 
practically over the entire frequency range, a mode cor- 
responding to almost free propagation in an effective ho- 
mogeneous medium, with Bragg gaps opening up about 
LoS/ca ~ 1.5 at the edges of the Brillouin zone and about 
LuS/ca ~ 3.5 at the center of the Brillouin zone; and we 
would have also resonance bands developing about the 
I ~ and I = 1 resonances of the individual spheres. 
From Fig. ^ we expect a single I — resonance band 
of Ai symmetry about ujS/ca ~ 2.2; and three I ~ 1 
bands about tuS/ca ~ 3.9, of which one should be non- 
degenerate (Ai symmetry) and one doubly degenerate 
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FIG. 7: (a) The phononic frequency band structure normal 
to the (111) surface of a fee crystal of silica spheres in air 
with / = 74%. The solid and dashed lines refer to bands of 
Ai and A3 symmetry, respectively, (b) Transmittance of an 
acoustic wave incident normally on a slab of the above crystal 
consisting of 32 planes of spheres parallel to the (111) surface. 
The positions of the first two resonant modes of a single silica 
sphere in air are indicated in the margin. 



(A3 symmetry). We expect a degree of hybridization be- 
tween bands of the same (Ai) symmetry, and this natu- 
rally leads to the bands shown in Fig. 0(a), except for the 
one flat band, a doubly degenerate band of A3 symme- 
try, about ijjS/ca = 3.426, which appears to have a very 
different physical origin (see below). 

In the long-wavelength limit {uj — > 0) we obtain a lin- 
ear dispersion curve, the slope of which gives an effective 
velocity of sound in the composite medium, c = 278 m/s, 
which is in very good agreement with result— c = 291 
m/s, of the effective-medium approximation. til We note 
that this agreement is greatly improved at lower values 
of the volume-filling fraction /. It is worth noting, also, 
that the effective velocity of sound in a system of sil- 
ica spheres in air is smaller than in air because of the 
high density contrast between silica and air. Indeed, in 
the limit pa << Ps, the effective- medium approximation 
gives practically identical results with the simple expres- 
sion c/cq = [2/ (2-f /)]^/^, which gives c < Cq at any value 
of/. 

It seems that the modes corresponding to the fiat 
band about ujS/ca — 3.426, which are deaf modes, are 
those of a wavefield highly concentrated between consec- 
utive planes of spheres, with very little interaction be- 
tween neighbor regions of high concentration. In order 
to demonstrate the above, we looked for the eigenmodes 
of the acoustic field, for k|| = 0, in a slab of TV = 2, 4, 8, 16 
planes of spheres, in the manner described for tiie corre- 
sponding problem of the electromagnetic field,Ea over a 
frequency range about the said fiat band of Fig. 0(a). Be- 
cause of the two-dimensional translation symmetry of the 
slab, the wavefield, in the air regions, can be expanded 
[as in Eq. (m)] in a series of plane waves with wavevec- 
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FIG. 8: The deaf band about ujS/ca = 3.426 of Fig. §(a) 
(solid line). The doubly degenerate eigenfrequencies of the 
corresponding bound states in slabs of A'' = 2, 4, 8, 16 planes 
of spheres are noted by circles. 



where it couples with an incident wave of the same ky 
leading to measurable transmittance. 

In Fig. 1^, we present the projection of the frequency 
band structure of the acoustic field of the phononic 
crystal under consideration, on the symmetry lines of 




r KM r 



tors Kg ~ ^g, ± \{uj/caY — g] "^^^^ ; now, as it turns out, 
for the eigenmodes in question the term corresponding 
to g = in the above plane- wave expansion vanishes, 
which implies that the wavefield decreases exponentially 
on both sides out of the slab, since ui/ca < |g| for g 7^ 0. 
We can say that these eigenmodes correspond to bound 
(in the z direction) states of the acoustic field. Our re- 
sults are shown in Fig. ^. For a single pair of planes 
{N = 2) we obtain a pair of such modes correspond- 
ing to a doubly degenerate eigenfrequency. For three 
pairs of planes {N = 4) we obtain three pairs of doubly 
degenerate eigenfrequencies, and similarly in each case 
we obtain a pair of doubly degenerate eigenfrequencies 
for every pair of planes in the slab. The eigenfrequen- 
cies corresponding to the bound states of the different 
slabs in Fig. ^ have been plotted against values of the 
reduced wavenumber k^i = iri/Nd, i — 1,2, . . . N — 1, for 
the slab consisting of N planes, to show how this band 
of states develops with increasing N. In every case, d 
equals the separation of two consecutive fee (111) planes 
in the infinite crystal, so that Nd equals the thickness of 
the slab. We note that the eigenfrequencies of the finite 
slabs coincide with the corresponding eigenfrequencies on 
the dispersion curves (solid lines in Fig. H) of the infinite 
crystal. The question now arises as to whether one can 
see experimentally this band. The corresponding modes 
do not couple to an incident wave. This is demonstrated 
quite clearly in Fig. |7|(b), which shows the transmission 
coefficient of a wave incident normally on a slab of the 
material consisting of 32 planes of spheres. We can see 
that the transmission coefficient vanishes at and about 
ojS/ca — 3.426 where the band under consideration ex- 
ists. However, we have verified that the said band sur- 
vives for k|| ^ (at least in the neighborhood of k|| — 0) 



FIG. 9: Projection of the phononic frequency band structure 
of a fee crystal of silica spheres in air, with / = 74%, on the 
SBZ of the fee (111) surface, along the symmetry lines shown 
in the inset. Propagating waves in the air about a slab of 
the crystal exist for frequencies above a threshold value (a 
function of ky) ujinf = Ca |k||| denoted by the dashed line. 



the SBZ of the fee (111) surface. The shaded regions 
extend over the frequency bands of the acoustic field: 
at any one frequency within a shaded region, for a 
given k||, there exists at least one propagating acoustic 
mode in the infinite crystal. The blank areas correspond 
to frequency gaps. We note that knowing the modes 
with k|| in the shaded area (FKM) of the SBZ and 
—Tr/d<kz < n/d is sufficient for a complete description 
of all the modes in the infinite crystal. The modes in the 
remaining of the reduced k space are obtained through 
symmetry. One clearly sees that for the opal under 
consideration, and unlike the colloidal crystal studied 
in Section ||, one obtains an omnidirectional frequency 
gap extending from ujS/ca = 1-59 to toS/ca — 1-95. We 
verified that this is indeed so by calculating the band 
structure at a sufficient number of ky points in the SBZ. 
Finally, we should note the existence of a very narrow 
omnidirectional gap at higher frequencies, extending 
from LuS/ca = 3.30 to LoS/ca = 3.34. 
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